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One of the key requirements for good thermoelectric materials is a low lattice thermal conductivity 
Here we present a combined neutron scattering and theoretical investigation of the lattice dynamics 
in the type I clathrate system Ba-Ge-Ni, which fulfills this requirement. We observe a strong 
hybridization between phonons of the Ba guest atoms and acoustic phonons of the Ge-Ni host 
structure over a wide region of the Brillouin zone which is in contrast with the frequently adopted 
picture of isolated Ba atoms in Ge-Ni host cages. It occurs without a strong decrease of the 
acoustic phonon lifetime which contradicts the usual assumption of strong anharmonic phonon- 
phonon scattering processes. Within the framework of ab-intio density functional theory calculations 
we interpret these hybridizations as a series of anti-crossings which act as a low pass filter, preventing 
the propagation of acoustic phonons. To highlight the effect of such a phononic low pass filter on the 
thermal transport, we compute the contribution of acoustic phonons to the thermal conductivity of 
BasGe4oNi6 and compare it to those of pure Ge and a Ge46 empty-cage model system. 



A central issue in thermoelectrics research is to find 
materials that generate a high electromotive force un- 
der a small applied thermal gradient. The best thermo- 
electric materials (TE) should therefore combine a low 
lattice thermal conductivity with high electrical conduc- 
tivity and thermopower. An efficient and economic way 
to reduce the thermal conductivity of bulk TE materials, 
without degrading their electronic properties, is to take 
advantage of inelastic resonant scattering between heat- 
carrying acoustic and non-propagative phonons, arising 
from isolated impurities with internal oscillator degrees 
of freedom. The librational vibration of molecules in- 
corporated in crystal structures is a well-known example 
of such scattering centers and was formerly observed in 
KC1 ionic crystals, doped with anionic molecules pQ, or 
in the filled water nanocages of clathrate hydrates |2j [3] . 
In TE clathrates and skutterudites, in which the phonon 
dispersions were recently investigated by means of in- 
elastic neutrons scattering (INS) [4H6| [TU]. the localized 
resonators are so-called rattling phonon modes of loosely 
bonded guest atoms in oversized atomic cages. The effect 
of inelastic resonant scattering [TTJ [12] was described the- 
oretically, using the relaxation-time approximation of the 
Boltzmann equation, in which the lattice thermal con- 
ductivity, ft/,, for a cubic crystal can be expressed as 

PUJmax 

kl = 1/3 / C v (cu)v(uj) 2 r(uj)po(uj)duj . (1) 
Jo 

Here C v (uo) is the phonon specific heat per unit volume 
and unit frequency, v(uj) is the mean group velocity for 



phonon frequency lj, t(uu) is the mean time between 
collisions that destroy the heat current for all phonons of 
frequency lj and po(u) is the normalized vibrational den- 
sity of states. The most common approach to evaluate 
the scattering of acoustic phonons by a resonator with 
frequency ljq is to add a phenomenological relaxation 
rate, t^ 1 (cj) = C ^jrz^p , where C is a constant propor- 
tional to the concentration of the resonant defects and 
to the strength of the resonator-lattice coupling. This 
has been used to fit experimental data for clathrates 
and skutterudites p~3j [14] . At a microscopic level, this 
implies the existence of a coupling and an anti-crossing 
behavior between the non-dispersive phonons of the 
isolated resonators and the acoustic phonon branches 
leading to the opening of a gap surrounding uo in the 
dispersion of the acoustic branch [17] . INS measurements 
on the clathrate BagGesoGaie [SJE] and the skutterudite 
CeRu4Sbi2 [4] evidence such an anti-crossing. However, 
in these previous studies, the authors conclude that 
the remarkably low k values of these materials cannot 
fully be ascribed to the effect of guest vibrations on 
lifetime t(lj) or group velocity v(uj) of the acoustic 
phonon modes (see also [H]). On the other hand, the 
picture of isolated guest atoms in a host cage and the 
microscopic mechanism responsible for the scattering 
of heat carriers were revisited by INS studies of the 
phonon dynamics in polycrystals of the skutterudites 
La/CeFe4Sbi2 [8 ]. The authors demonstrated that the 
guest atoms are coherently coupled with the host-lattice 
dynamics and associated with low energy optical phonon 
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modes which are characterized by a small group velocity. 
It was then proposed that these low frequency phonon 
branches open new paths for the occurrence of Umklapp 
scattering of heat-carrying acoustic phonons, that is, 
the backscattering of heat carriers by an elastic or 
a multiphonon process mediated by a lattice vector 
[4j El [8] . This mechanism obeys momentum and energy 
preservation; hence, it is in contrast to inelastic resonant 
scattering associated with energy dissipation at isolated 
atoms. The reason to invoke such ph-ph processes in 
complex unit cell materials is the presence of flat optical 
phonon modes at low energies which are expected to 
enhance ph-ph scattering processes. However, to date 
experimental evidence for this claim is missing, - it 
would require measuring the corresponding broadening 
of the energy width of the phonons involved. Yet, the 
available INS measurements on single crystals [4j |6] as 
well as our data demonstrate that the finite lifetime of 
acoustic phonons is not small enough to explain such a 
low ft, meaning that no fingerprints of strong Umklapp 
scattering could be evidenced. Thus, the microscopic 
origin for the low n remains elusive. In particular the 
respective roles played by guest atoms and unit cell 
complexity are ambiguous [T8] . 

In this paper, we present in section I the experimental 
mapping of the transverse acoustic (TA) and the guest 
phonon dispersions in a single crystal of the clathrate 
Ba8Ge42.1Ni3.5D0. 4 [24] by INS, covering the entire first 
Brillouin zone (BZ) and a wide energy range. The high 
quality of the data allows us to analyse the relative 
change in intensity of acoustic and optical phonons. We 
evidence a gapless dispersion of the TA phonons in ab- 
sence of a strong broadening of their energy widths when 
approaching the energy of the optical phonons. Moreover 
the existence of a mechanism of spectral weight trans- 
fer between acoustic and optical phonons resulting from 
strong hybridizations is visible over a wide region in recip- 
rocal space. In section II, we interpret our experimental 
results in the framework of ab-intio DFT calculations 
comparing the lattice dynamics of pure Ge, the empty 
Ge46 framework structure and the Ba8Ge4oNi6 clathrate. 
By a direct comparison with the experimental findings, 
we highlight the roles played by unit cell complexity and 
guest atoms for the lattice dynamics, respectively. We 
conclude that the low lattice thermal conductivity is due 
to a phononic filter effect. 

I. INELASTIC NEUTRON SCATTERING 

Subsection a) of this paragraph, contains details on 
the investigated sample, while in subsections b) and c) 
we provide the basic theoretical background of phonon 
measurements by INS and explain how our study differs 
from previous measurements done on similar materials. 



In subsection d) we report our results together with our 
analysis and interpretation of the data. Finally, in sub- 
section e) we discuss technical details with respect to the 
validity and reliability of our data analysis. 

a) Sample preparation and physical properties 

A single crystal of the clathrate Ba8Ni3.5Ge42.iDo.4 
(space group Pm3n, a = 10.798(2) A, 0= 8 mm and 
height = 30 mm), as shown in Fig. [I] a), was grown from 
the melt, using the Bridgman technique. A detailed 
analysis of its crystal structure and its thermoelectric 
properties were reported in [24]. This study revealed 
Ba8Ni3.5Ge42.1nc4 to be a thermoelectric (TE) n-type 
metal with a relevant value for the dimensionless TE fig- 
ure of merit ZT at high temperatures, thanks to its very 
low thermal conductivity of about k ~ 1.2 W/mK at 
300K. 

b) Neutron intensity, phonon dynamical structure 
factor 

In case of coherent inelastic nuclear scattering by a 
phonon of branch j, with energy uj q j and polarization 
vector £?, , the neutron diffusion function is written as: 

S ph (Q,uj) = n(Lj) ^ Q ^ 2 S(lj - uj qJ )S(Q - q- G) (2) 

U q,3 

where n(u) = i- exp (- hoj /k B T) is tne Bose factor > Q = Q+ 
G is the scattering vector given by the nearest reciprocal 
lattice vector G and the phonon wave vector q. Finally 
F 3 D (Q) is the dynamical structure factor (DSF) defined 
as 

f d(Q) = ^e- w ^^e i ^{Q.C J (Q)} (3) 

where 6 i? f^, Mi and Wi{Q) are coherent scattering 
length, fractional coordinates, mass and Debye- Waller 
factor of the i'th element, respectively. The scalar prod- 
uct {Q-Cj(Q)} contains the phonon polarization and can 
thus be used to distinguish longitudinal and transversal 
phonon modes by choosing the the appropriate combina- 
tion of phonon wave vector q and reciprocal lattice vector 
G. The polarization vectors of longitudinal and transver- 
sal phonons are parallel and perpendicular to the phonon 
wave vector g, respectively. In our case, as sketched 
in Fig.JT] c), measurements were conducted around the 
G = (006) reflection, scanning along gjno]- In this set- 
ting longitudinal phonons, polarized parallel to <f, are 
perpendicular to G and therefore only transverse compo- 
nents are selected. 

It is then worth to denote the diffusion function in 
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FIG. 1: (Color online) a): The measured single crystal 
of Ba8Ni3.5Ge42.1D0. 4 [24] b): Elastic neutron scans around 
the Bragg peak (600) along the transverse, [110], and lon- 
gitudinal, [100], directions, c): Bragg peak intensities in 
the {[100]; [011]} plane, calculated for a neutron wave vector 
k = 2.662A -1 . The pointsize is proportional to the nuclear 
elastic structure factor. The intensity of the strongest Bragg 
peak (600) is then compared to those used in [5] and [6]. As 
sketched and explained in subsection Lb), INS measurements 
of TA phonons are obtained by scanning the reduced wave 
vector, q, along the [011] direction from the Bragg peak (600). 

the long wavelength limit, for \q\ « \Q\. It can be 
shown that for acoustic phonon modes, close to the Bril- 
louin zone center, all atoms are vibrating in phase and 
with equal displacements £j(Q) / W^i- Consequently, for 
T » cj, the diffusion function, Eq. [2j becomes propor- 
tional to the nuclear structure factor (i.e. the Bragg in- 
tensity) : 

SMQ^)*G 2 ^4^ (4) 

Thus, the intensity is proportional to the nuclear struc- 
ture factor of the respective Bragg peak. Moreover, close 
to the Brillouin zone center and for pure acoustic phonons 
S p h{Q,u) x Uj(q) is constant. 

c) Measurement strategy and data analysis 

INS experiments were carried out on the thermal triple 
axis spectrometer 2T at the Orphee reactor of the Labo- 
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FIG. 2: (Color online) a): Energy scans close to the Bragg 
peaks Q=(600) or Q=(222), at wave vector qy = 0.26lA _1 
in transverse geometry, b): Energy scans measured at the 
constant wave vectors Q = (611), open circles, and Q = 
(5, 2.1, 2.1), open squares. Three Gaussian profiles (blue lines) 
centered at the energies labeled by E±, E2 and E3 were used 
to fit the total phonons line shape (black line). The instru- 
mental energy resolution, AEi n3 t=l.l meV, is indicated by 
the red line 



ratoire Leon Brillouin (LLB, France). All measurements 
were conducted at room temperature in fixed Q mode 
with a final wave vector of 2.662 A -1 . To eliminate 
higher-order harmonics graphite filters were used. The 
sample was mounted in a vacuum chamber and aligned in 
the plane {[100]; [011]} such that wave vectors of the form 
Q=27r/a(hkk) were accessible as depicted in Fig. [l]c. 
Wide elastics scans around the Bragg peak (600) are 
shown in the Fig. [T] The bulk mosaicity is of 0.75°. 
In the long wavelenght limit the acoustic phonon inten- 
sity is directly proportional to the Bragg intensity as 
outlined above, see Eq. [4] Therefore acoustic phonons 
have to be measured near strong nuclear Bragg reflec- 
tions. The main technical difference between the present 
and previous INS experiments on clathrates is that we 
could access a higher Q-range and thus measurements 
around the strongest Bragg peaks could be conducted. 
As sketched in Fig. [l]b, our measurement was conducted 
around the strongest reflection G = (600), following the 
transverse [Okk] direction, while previous INS studies of 
transverse phonon branches were carried out around less 
intense Bragg peaks (400), (222) and (330) [6]. As 
a consequence of the INS cross section detailed above, 
the TA phonons in our data appear much stronger than 
the guest phonon modes and with a significantly en- 
hanced resolution, when compared to the results in Refs. 
[5] and [6]. To emphasize this aspect, we compare in 
Fig. [2] a) transverse phonon spectra measured at the 
same wave vector |<fi|(0.26lA -1 ) either around the Bragg 
peaks Q=(600) or Q=(222). Another consequence of the 
measurement around the Bragg peak G = (600), which 
follows directly from energy and momentum conservation 
laws, is that our measurements go up to 14 meV while 
data in Ref. [6] stop at 7 meV. This is an important ad- 
vantage which allows us to follow the spectral weight dis- 
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tribution and transfer among the different acoustic and 
optic phonon branches. 

Thanks to the high data quality, we are then also able 
to perform an advanced and reliable analysis in order to 
extract properly the intrinsic energy width, the energy 
position and the dynamical structure factor, or norm, of 
the contributing phonons. In fact, realistic INS measure- 
ments on a conventional (mono-detector) neutron triple 
axis spectrometer mean scanning the (Q, oS) space step by 
step with a resolution function. The neutron intensity at 
a given point (Qo^o) is given by the convolution prod- 
uct of this instrumental resolution function, R(Q, cj), and 
the diffusion function, S(Q,uj), which contains all the in- 
formation of the dynamics in the sample probed: 

HQo^o) = S(Q, uj) (8) R(Q -Q ,uj- uj ) (5) 

The resolution function is a 4-dimensional function 
which couples the energy and the reciprocal space. A 
cut of this resolution function in the {qokk,w} plane is 
shown in Fig. [5] a). As one can see, the size of this 
function is not negligible with regard to the phonon 
branch dispersion. Therefore, in order to extract 
properly the intrinsic position, width and DSF of a 
phonon from the measured scattering profile, we use a 
fitting procedure which numerically convolutes a model 
phonon cross section with the spectrometer resolution. 
It enables us to take into account the local curvature of 
the phonon branch over the range of the instrumental 
resolution, i.e. the fact that the energy resolution for an 
acoustic phonon mode is not independent of its slope. 
Only this method allows for a reliable extraction of the 
characteristics of phonons from INS data. The DSF 
has so far not been discussed in any neutron work on 
clathrates. (For a similar analysis on the skutterudite 
CeRu4Sbl3 we refer to g|). 



(d) INS Results 

The energy scans, shown in Fig. [2]b), were performed 
far from any acoustic phonon branch at the wave vectors 
(611) and (5,2.1,2.1) to localize the low lying optical 
phonon modes. The energy spectra at (611) were ana- 
lyzed by using three Gaussian profiles centered at ener- 
gies of £i=5.5(0) meV, £ 2 =8.2(0) meV and £ 3 =10.8(5) 
meV with full width at half maximum (FWHM) val- 
ues of AE 1 ~ 2 meV, AE 2 ~ 2 meV and AE 3 ~ 3 
meV. Since the linewidth of these peaks is larger than 
the energy resolution, several phonon branches can be 
contained in a peak. Later, the peaks E\ and E 2 will be 
identified as stemming from localized Ba motions in the 
large Ba(Ge,Ni)24 cages, while the peak E% contains both 
optical phonon modes from the (Ge,Ni)46 host framework 
and guest phonon modes, in agreement with previous INS 
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FIG. 3: (Color online) Energy scans with constant Q = 
(600) + q(Ofcfc) in the first BZ, k <0.5. The solid black lines 
are fits as explained in the text. Red lines and filled areas 
depict contributions of TA phonons while blue lines are rep- 
resentatives of the peaks at E\ , E 2 and E 3 . For scans labeled 
by -E the negative energy region is shown (Stokes side). 



[EIHHEHIEI], Raman [28] and theoretical [1[25] reports. 
In a next step the TA phonon branch was investigated 
by performing energy scans around the Bragg peak (600) 
for a set of reduced phonon wave vectors, q = (Ohk), 
covering an entire BZ (0 < k < 1). Raw data are shown 
in Fig. [3] for k <0.5 and in Fig. [4] for k >0.5. The model 
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FIG. 4: (Color online) Energy scans with constant Q = 
(600) + q(Ofcfc) for k > 0.5. The solid black lines are fits as 
explained in the text. Red lines and filled areas depict con- 
tributions of TA phonons while blue lines are representatives 
of the peaks at E\ , E2 and E3 . 



phonon cross section, applied for our fitting procedure, 
consists of a damped harmonic oscillator for TA phonons 
and three Gaussian functions (centered at energies Ei, 
E2 and £^3) for the low- lying optical phonons as before- 
hand identified. In Fig. [2] red lines and shaded areas 
depict the contributions of the TA phonons, while the 
blue lines represent the three optic-like excitations at E\ , 
E 2 and E%. The deconvoluted energy positions and dy- 
namical structure factors (DSF) of both the TA phonons 
and the three optic-like excitations are shown in Fig. [5] as 
functions of q. 

For k <0.5, the phonon spectra consist of a main peak, 
stemming from the TA phonon branch (red lines in 
Fig. [3]), which appears at energies distinctly lower than 
the optical excitations. As shown in Fig. [5] a), the energy 
positions of these peaks describe a dispersive branch be- 
low the optical continuum. This branch does not cross 
any optical phonon mode and no gap is observed follow- 
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FIG. 5: (Color online) a) Energy dispersions of the TA and 
the three optic- like (Ei, E2 and E3) phonons. The size of 
the error bars represent the phonon linewidth (FWHM). The 
ellipse is a cut of the instrumental resolution volume (axis 
sizes: 0.03 A" 1 (0.034 rlu) along the (Okk) axis and of 0.24 
meV along the energy axis), b) DSF as function of q in the 
direction (Okk). 



ing its dispersion. For wave vectors below about &=0.25, 
the energy width of the TA modes is resolution limited 
and their DSF fluctuates around a constant value. In the 
limit uj — » 0, the slope of the TA phonon branch yields 
the transverse sound velocity z^ 011 ^=2530 m/s, compa- 
rable to that one of BagGesoGaie [6]. The dispersion of 
the TA mode deviates from linearity for energies above 
about 4 meV (~ 46 K). For 0.25 < k < 0.4, although 
there are still well defined TA phonons with resolution 
limited energy widths, one observes an abrupt decrease 
of their DSF. The DSF of TA phonons drops by more 
than 80% between k = 0.25 and k = 0.4, as seen in Fig. [5] 
b). The analysis of the experimental data becomes more 
delicate for wave vectors above k=0A where intensities 
and energy positions of the TA branch and the low lying 
optical excitations come closer. However, as can be ob- 
served in the raw data in Fig. [3] as well as in Fig. [6] where 
only the phonon spectrum measured at k=0.4062 is de- 
picted, there are always well defined TA phonon peaks 
located at energies below the optical excitations. The 
profile of these peaks is continuous until a break at about 
the energy E x of the low lying optical excitation. While 
the TA phonon branch fades away, the intensities of the 
guest phonon peaks at E\ and E 2 change (see Fig.|5]b)). 
The DSF of the guest peak at E\ increases from a con- 
stant value for k <0.25, reaches its maximum at g=0.43 
A -1 (fc=0.5) and then returns to roughly the same con- 
stant value. 

A similar ^-dependence of the DSF is observed for the 
guest peak at E 2 with the maximum slightly shifted to 
a higher q of about 0.54 A -1 (fc=0.62). For the peak 
at Es no sizable change is observed over the entire q- 
range. Note that the qualitative changes of the guest 
DSFs are clearly visible on the raw data in Figs. [3] and 
[4] (compare e.g. data measured at /c=0.34 to those ob- 
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tained at /c=0.46, k=0.62 and fc=0.87). As evidenced in 
Fig.[5]b), the ^-dependence of the DSF seems to arise 
from a transfer of spectral weight between TA phonons 
and guest vibrations. As k increases, the spectral weight 
of the TA phonons is first transferred to the first guest 
phonon peak at E\ and then to the second one at E2. 
Interestingly, similar effects have been observed in qua- 
sicrystals and large unit cell crystals [191 [20] . The energy 
of the guest phonons at E\ increases significantly during 
this process, shifting from 5.4 meV at the zone center to 
6 meV at /c=0.5 without any sizable change in its energy 
width (see Fig. [5] a)). 

To conclude, firstly, our data put in evidence a continu- 
ous dispersion of the TA phonon branch throughout the 
whole BZ zone with a rapid bending together with a dras- 
tic decrease of the spectral weight as k increases beyond 
0.25 or for energies higher than 4 meV, as summarized 
in Fig. [5] a) and b). The most surprising and relevant ef- 
fect observed is precisely that the spectral weight of TA 
phonons decreases abruptly without a strong broaden- 
ing of their energy profiles. In others words, within our 
resolution, no finite energy damping or finite lifetime is 
needed to fit the TA phonon profiles. Note that the lack 
of line broadening in the deconvoluted data is not meant 
to demonstrate the complete lack of anharmonic phonon 
scattering processes, it only emphasizes that anharmonic 
effects are by far not as big as assumed in previous stud- 
ies. Therefore it becomes evident that it is not correct to 
consider a strong decrease of the relaxation time r(u) as 
function of energy in Eq. [I] to account for the effect of 
guest phonon modes. Moreover, the TA phonon branch 
disperses continously without any fingerprint of a gap 
induced by a crossing with the low lying optical phonon 
branch in contrast to what is reported by [6] . This makes 
it obvious that adding resonant like relaxation times in 
Eq.[l]to account for the effect of vibrations of the guest 
atoms is incorrect. 

Secondly, our data highlight a mechanism of spectral 
weight transfer between TA and guest phonon modes 
which occurs over a wide region of the BZ, from k=0.25 
to k=0.5. Such transfer of spectral weight was observed 
in the skutterudite CeRu4Sbi2 @|. Note that along a 
given direction across the BZ, neither an intensity change 
nor an energy dispersion would be expected for guest 
phonons, if they were assumed to be isolated oscillators. 
Since we observe changes of the spectral weight of both 
Ei and E2 guest phonons as well as a dispersion of the 
peak Ei, the picture of a "freely rattling" guest atom in 
a host cage is not applicable in the case of clathrates, as 
previously concluded from INS experiments on skutteru- 
dites 0. 
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FIG. 6: (Color online) a) Fit to the scan at /c=0.4062 rlu, 
obtained by forcing a reduction of the intensity of the optic 
phonon mode at E\ (shaded in blue). The fits provide the 
upper limit of the intensity of the acoustic phonon as indicated 
by error bars in Fig. [5] a), b) green area represents the 
simulated profile of the phonon measured at /c=0.3749 rlu 
with a finite energy width of Y — 1.2 meV. 



(e) Validity and reliability of the data analysis 

Crucial for our data analysis is the ability to differen- 
tiate between the acoustic and optical parts of the mea- 
sured spectra even in the region where both contributions 
are superimposed, i.e. for k >0.4. The TA phonons have 
lost already 80 % of spectral weigth at k — 0.4, see Fig. [5] 
b), with no evidence of any broadening in energy. Thus 
the main effect of spectral weigth transfer between the 
TA and optical phonons occurs then well below k =0.4 
where the precision of the fits is excellent. We now dis- 
cuss the data analysis of the energy scan measured at 
/c=0.4062, shown in Fig. [6j in which the distance in en- 
ergy between the TA phonon and the optical phonon at 
Ei is less than 1 meV. In our data analysis, the surpris- 
ing result is that, within our resolution, no finite energy 
width is needed to fit the profile of the TA phonons. The 
observed energy width arises simply from the integration 
of all the phonon modes at 0.4062 ± 5k contained in the 
instrumental resolution function. However, one can try 
to force the occurrence of a finite energy width by forcing 
the intensity of the E\ peak contribution to decrease as 
compared to the fitting model proposed in Fig. [3] The 
result, shown in Fig. [6] a), clearly yields a poorer fit with 
an intrinsic energy width of 0.45 meV and an increase 
of about 25% of the spectral weigth of the TA mode. 
The resulting limits for the TA phonon intensity are in- 
dicated by the red line on the y-axis in Fig.[5]b). We thus 
clearly evidence that, independtly of the fitting model, 
the acoustic phonon looses abruptly its intensity when 
approaching E\. 

As demonstrated above, the upper limit for the TA 
phonon energy width is roughly 0.45 meV. In a sim- 
ple Debye approach of the thermal conductivity based 
on Eq. [l] and on our measured TA phonons disper- 



7 




k (rlu) k (rlu) k (rlu) 



FIG. 7: (Color online) Calculated dispersion curves along the 
reciprocal space direction (6kk) for a) BasGe4oNi6, b) Ge46 
and c) diamond Ge. Modes are color coded with respect to 
their PR. While green means a high PR, red and black stand 
for intermediate and low PR, respectively (see text). The 
Energy axis of panel a) is rescaled by a factor 1.4 to match 
the experimental data. 

sion (see below), we estimate that the relaxation time 
that would be required to reproduce the low experimen- 
tal value of the thermal conductivity in our material is 
of about r ~ 0.2 ps corresponding to a phonon energy 
width of about 3.3 meV. In Fig.[6]b), we simulate the 
profile of the TA phonon measured at &=0.3749 with a 
finite energy width of 1.2 meV, dashed line, and we com- 
pare it with the result of our best fit. It evidences that 
the main microscopic mechanism responsible for the low 
thermal conductivity cannot be a large relaxation time. 

II. DFT CALCULATIONS 

To deepen our understanding of the INS measurements 
and of the coupling mechanism between the guest and 
host phonons, we have performed an ab initio density 
functional theory (DFT) study of the lattice dynamics 
and the INS cross section in the cubic high-symmetry 
configuration BagGe^Nie and an empty Ge46 framework 
model system. 

Ab-initio density functional theory (DFT) calculations 
were conducted using the DFT code VASP [29] apply- 
ing the projector augmented wave generalized gradient 
method (PAW-GGA). First the unit cell of Ba 8 Ge 40 Ni 6 
was relaxed to its electronic ground state, using high pre- 
cision settings with a 4x4x4 k-point mesh and a con- 
vergence criterion of residual forces of less than 10 -4 
eV/A. After reaching the ground state, symmetry non- 
equivalent displacements were introduced and the re- 
sulting Hellmann-Feynman forces were calculated and 
used to determine the dynamical matrix by use of the 
PHONON package [30]. Different displacements, ranging 
from 0.03 A to 0.08 A, were investigated, yet no signif- 
icant differences in the phonon properties could be ob- 
served. The presented calculations were then conducted 
for displacements of 0.05 A. By solving the dynamical 
matrix for selected q- values, phonon eigenfrequencies and 
corresponding eigenvectors can be determined within the 
harmonic approximation. From these eigenvectors and 




k (rlu) k (rlu) k (rlu) 

FIG. 8: (Color online) Calculated dispersion curves as in Fig. 
[7] over the full energy range. 

eigenvalues it is then possible to calculate diffusion func- 
tion and DSF following Eqns. Q and 

a) Dispersion curves 

Figures]?] and [8] depict a comparison of the phonon dis- 
persion curves of BagGe^Nie and Ge46 along the recip- 
rocal space direction [Okk]. As for other complex systems 
[22] , a rescaling of the energy axis is necessary to achieve 
good agreement between the absolute energy values of 
calculated and experimentally determined acoustic and 
low- lying optical branches in BagGe^Nie. The rescaling 
parameter of 1.4 is relatively large, pointing to an under- 
estimated Ba-Ge interaction within the DFT approach. 
Colors in Figs. [7] and [8] represent the phonon participation 
ratio (PR), defined for an eigenmode j as [23 

P»= E^mH l N T.^Mf~ ■ ( 6 ) 

\i=l 1 / i=l 1 

The PR is calculated from the phonon eigenvector £j (q) 
of a given mode and the masses of the respective atoms. 
A PR close to 1 (green) means coherent displacements 
for all atoms of the same species and thus a collective 
excitation. Localized excitations, such as rattling modes, 
have a PR close to (black), meaning that only a few 
atoms have large displacements, while the others are 
only slightly or not at all displaced. 

The comparison of the two cage structures with the 
diamond Ge enables to disentengle different reasons 
for a reduction of k>i. If we first compare Ge46 and 
diamond Ge we find that many phonon modes in the 
complex structure are less dispersive and evidence a low 
PR. Interestingly this corresponds nicely to the results 
obtained in reference [16 , where already the clathrate 
framework is found to reduce ^. While neither pure Ge 
nore the empty Ge46 structure evidence modes with a 
very low PR, the flat dispersionless branches between 
6 and 8 meV in Ba8Ge4oNie clearly do and can thus 
be identified as Ba rattling modes. These branches 
have strong impact on the acoustic phonon modes. 
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First, they force the acoustic branch to bend over at 
lower energies, thus decreasing the group velocity of 
the acoustic phonons. Furthermore, they successively 
hybridize with the acoustic branch and lead to an 
important reduction of the PR of the hybridized acoustic 
waves, which will contribute to the dissipation of the 
heat current through the material. More than a simple 
ant i- crossing, the rattling modes act as a low pass filter 
for propagating phonons which prevents propagation of 
heat carrying acoustic phonons with energies higher than 
the Ba rattling energies. Filtered phonons are not only 
dispersionless modes with small group velocities but are 
also characterized by a low or intermediate PR, meaning 
that only few atoms are actually involved. Unfiltered 
phonons are unaffected by the rattling phonons and 
their efficiency to propagate heat depends on the usual 
anharmonic scattering processes. 

If we now moreover investigate the simulated INS scat- 
tering profiles, a sequence of spectral weight transfer from 
acoustic phonons to the different rattling modes is visible, 
as was found for the experimental data of BagGe^Nig, 
as k is increased (see Fig. |9|. However, these effects are 
clearly more pronounced in the experiment, in which hy- 
bridization starts at lower k values and leads to a faster 
decrease of the DSF of the TA phonon modes. This might 
be attributed to the fact that, in our simulation, we use 
a model which is free of disorder. Indeed it can be shown 
that the Ni content strongly influences the cage size and 
the related Ba rattling frequencies. 
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FIG. 9: Constant q scans in the first BZ, determined for 
the same q values as measured experimentally. The finite 
linewidths is obtained by convolution with a Gaussian func- 
tion of 0.1 meV width. 



Fig. 10 is splitted in two branches on the bottom panel, 
thus confirming that disorder yields a broadening of the 
rattling modes. 



b) Disorder 



c) Phononic filter effect 



Due to the size limit which is imposed to ab-initio cal- 
culation, an investigation of disorder is rather difficult. 
In the Ba-Ge-Ni system, it is yet possible to create ide- 
alized structures with different Ni content, which still 
evidence high symmetry and can be described in a unit 
cell approach, thus enabling us to get some idea about 
the influence of structural disorder. The structures of 
BagGe^Nie and BasGe42Ni4 differ on the 6c position, 
which is the only site that exhibits disorder in the Ba- 
Ge-Ni clathrate structure. While it is fully occupied by 
Ni atoms in the case of Ba8Ge4oNi6, it contains four Ni 
and two Ge, which are symmetrically distributed in case 
of our Ba8Ge S 42Ni4 model system. Since the 6c posi- 
tion is located in the large 24-atom trapezohedron, this 
cage indeed is influenced by the Ni content. An increas- 
ing Ge and vacancy content on this position thus intro- 
duces distortions and breaks the symmetry of this cage. 
This symmetry-breaking is expected to broaden the fre- 
quency distribution of the encaged Ba rattling modes. 
And indeed by comparing the DSF of Ba8Ge4oNie and 



To highlight the effect of such a phononic low pass filter 
on the thermal transport, we compute the contribution of 
LA phonons, propagating along the [Oil] direction, to the 
thermal conductivity of pure Ge, Ge46 and Ba8Ge4oNi6- 
From our ab-initio calculations for these three systems, 
we extract po(oo) and v(u) for the chosen phonon branch. 
We then evaluate Eq. ([I]) by assuming a constant phonon 
lifetime, r{uo) = 12 ps, estimated from the upper limit of 
the experimentally evidenced linewidth of TA phonons 
in the Ba-Ge-Ni clathrate. Moreover we replace uOmax by 
a cutoff frequency Co> cn t, which is chosen to be the fre- 
quency at which the LA phonon branch looses its acous- 
tic character, i.e. the frequency at which the PR of the 
acoustic mode strongly decreases (change from green to 
red in Figs. [7] and [8|. The introduction of this cutoff fre- 
quency takes the effect of the rattling modes into accout, 
which are responsible for a further flattening of the opti- 
cal branches, as was discussed in subsection]] Figure 11 



shows the frequency dependence of ftj?^ 1 ' (cj) along a LA 
branch at room temperature, as well as the resulting k(T) 
in the inset. The reduced thermal conductivity of Ge46 
with respect to pure Ge shows that the cage structure 



BasGe42Ni4 this is evidenced (see Fig. 10). The low ly- 
ing Ba rattling mode which is seen on the top panel of 
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FIG. 10: DSF of a) Ba 8 Ge4oNi 6 and b) Ba 8 Ge 4 2Ni 4 along 
the direction (6,£,£). 



itself is significantly decreasing the lattice thermal con- 
ductivity. An additional reduction of n in BagGe^Nie 
arises from the low energy Ba rattling phonons which 
lower the cutoff energy. Note that if for pure Ge we use 
the same cutoff frequency as for Ge46, this results, de- 
spite the differences in v(cj), in an similarly low thermal 
conductivity, as can be seen from the dashed curve in 
Fig. [II] Thus, to lowest order, the thermal transport by 
acoustic phonons in the Ba-Ge-Ni-clathrates corresponds 
to that of a simple Debye system with the sound velocity 
of pure Ge and a cutoff frequency corresponding to the 
energy of the first rattling mode of the guest atoms. 
It is of course true that to fully calculate the thermal 
properties of a material optical phonon modes also have 
to be taken into account, therefore our model calculation 
has to be understood as a first order approximation. This 
is justified, since the main contribution to heat transport 
can be attributed to acoustic phonons, while the optical 
modes above the introduced cutoff frequency can be neg- 
elected in a first order approach due to smaller lifetimes 
and group velocities. 



5 10 15 20 25 30 
E (meV) 

FIG. 11: (Color online) Calculated frequency dependence of 
the lattice thermal conductivity per unit frequency for pure 
Ge, Ge46 and BasGe^Nie at room temperature along a LA 
branch of the dispersion curve. The dashed blue line shows 
kla(w) for pure Ge with the same cutoff that is applied for 
Ge46- The inset shows the resulting k(T) for the three systems 
obtained from the average over TA and LA branches of the 
(Okk) direction. The curves for Ge (black) are scaled by a 
factor 0.2 (main panel) and 0.1 (inset) for better readability. 



III. CONCLUSION 

We have performed a high resolution INS study of the 
lattice dynamics of one of the prototypal thermoelectric 
clathrate systems, Ba-Ge-Ni. Our data exclude an in- 
terpretation in terms of an isolated oscillator, since co- 
herent modes of guest and host system are evidenced. 
Furthermore our results point out that Umklapp scat- 
tering is not able to account for the reduced kl in the 
Ba-Ge-Ni system. The evidenced phonon lifetimes are at 
least an order of magnitude bigger than what would be 
expected if such a low kl was to be reached by Umk- 
lapp scattering processes. Therefore a new mechanism 
was suggested to explain the low lattice thermal conduc- 
tivity in clathrates - the phononic low-pass filtering of 
acoustic phonon modes, having its origin in both loosely 
bound guest atoms and structural complexity of the host 
framework. We are confident that these new insights will 
promot further studies in clathrates and other guest-host 
materials. 
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